Directional Interactions between Membrane Proteins 
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While modern structural biology has provided us with a rich and diverse picture of membrane 
proteins, the biological function of membrane proteins is often influenced by the mechanical prop- 
erties of the surrounding lipid bilayer. Here we develop an analytic methodology connecting the 
hydrophobic shape of membrane proteins to the cooperative function of membrane proteins induced 
by bilayer- mediated elastic interactions. Application of this methodology to mechanosensitive chan- 
nels shows that, in addition to protein separation and bilayer material properties, the sign and 
strength of elastic interactions, and associated cooperative gating characteristics, can depend on the 
protein shape and orientation. Our approach predicts how elastic interactions affect the molecular 
structure, organization, and biological function of proteins in crowded membranes. 
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Cell membranes exhibit a complex organization of 
lipids and membrane proteins [l], Q and play an integral 
role in many cellular processes. The functional proper- 
ties of membrane proteins are not purely determined by 
protein structure but, rather, membrane proteins act in 
concert with the surrounding lipid bilayer [ll-Q- I n P ar_ 
ticular, the hydrophobic regions of membrane proteins 
couple to the hydrophobic regions of lipid bilayers [H- 
[ll| . The energetic cost of the resulting membrane defor- 
mations can be captured by an elastic model 12-18[ in 
which the lipid bilayer is described as an elastic medium 
and the membrane proteins are regarded as rigid mem- 
brane inclusions. Neighboring membrane proteins may 
induce overlapping deformation fields of the bilayer mem- 
brane, yielding long-range interactions between mem- 
brane proteins 16|, [l9|. Experiments have implicated 



such membrane-mediated interactions in the clustering 
of a number of different membrane proteins [l9l-[25|. 

Theoretical studies of membrane-mediated interac- 
tions between proteins [Tg|, [261-37] have mainly focused 
on the idealized case of cylindrical or conical membrane 
inclusions at larg e separations. However, cell membranes 
are crowded [38l440j , with the size and spacing of mem- 
brane proteins both being of the order of a few nanome- 
ters. In addition, structural studies [4l|, |42| have shown 
that membrane proteins exhibit a rich variety of hy- 
drophobic shapes. How are the shapes of neighboring 
membrane proteins reflected in the structure of elastic bi- 
layer deformations, and what are the resulting directional 
interactions between membrane proteins at the small sep- 
arations most relevant for cell membranes? In this Letter 
we describe an analytic methodology for estimating di- 
rectional interactions between membrane proteins for ar- 
bitrary protein separations. We first develop our general 
methodology, and then apply this methodology to the 
mechanosensitive channel of large conductance (MscL) 
43l-l47| . Our calculation of the interaction potentials and 
cooperative gating curves for the proposed pentameric 



structures of MscL 41, 47 -57} shows that, similarly as 



lipid bilayer material properties (12M18L |25|, l54H57j , elas- 
tic interactions can alter the structure and function of 
proteins in crowded membranes. 

Following the standard framework of membrane elas- 
ticity [58l46Q| , we describe the positions of the inner and 
outer lipid bilayer leaflets within the Monge representa- 
tion throu gh t he functions h+(x,y) and h-(x,y), which 
determine [61( the midplane and thickness deformation 
fields governing bilayer membranes. Since thickness de- 
formations are thought to be dominant for MscL gating 
62^651, we illustrate our methodology for the thickness 
deformation field u(x,y) = - [h+(x,y) — h-(x,y) — 2a], 
where 2a is the hydrophobic thickness of the unperturbed 
lipid bilayer, with the elastic energy 00,158^671 

G = \ J dxdy{K b (V 2 u) 2 + K t (^) 2 + T [2* + (V«) 2 ] } , 

(1) 

where is the bending rigidity of the lipid bilayer, K t is 
the stiffness associated with thickness deformations, and 
r is the membrane tension. Midplane deformations de- 
couple to leading order from thickness deformations [6l| , 
and are described [HsMgHI] by an energy functional similar 
to Eq. (pQ). For generality we allow for the two tension 
terms u/a and (Vu) in Eq. (pQ), which capture the effects 
of membrane tension on lipid surface area [go|, @, [fill and 
on membrane undulations [58l-[63|, respectively. 

The Euler-Lagrange equation associated with Eq. (pQ) 
is given by 

{V 2 -v+) (V 2 -v-)u = 0, (2) 
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where 2K h v± = r ± (r 2 - 4K b K t /a 2 ) and u = u + 
ra/K t . The solution of Eq. (j2j) for a sing le cylindrical 
membrane inclusion of radius is [12L 



(3) 



where ri and 6i are polar coordinates with the cen- 
ter of inclusion i as the origin, the Fourier-Bessel se- 
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ries f*(ri,0i) = A^ K (v^± r O + E^=i (A,n + S i>n ), 
in which A,n = A± n K n (y/v±ri) cos nOi and 23^ = 
Bf n K n (-^j7±ri) sin n^, if n are modified Bessel functions 
of the second kind, and we have assumed that mem- 
brane deformations decay away from a single membrane 
inclusion [66j. The coefficients Af n and Bf n are de- 
termined by the boundary conditions at the bilayer- 
inclusion interface for which, consistent with previous 
studies [12|, [26|, I62l - l67j , we use 



u{r h 6i)\ r .- R . = Ui(6i), 
n.V U (r,A)| n=fli = Ufa), 



(4) 
(5) 



where n is the unit normal vector along the bilayer- 
inclusion interface. 



Following Refs. [27|, [3l| we construct the solution of 
Eq. ([2]) appropriate for two cylindrical membrane inclu- 
sions of radius R\ and R2 using the ansatz 



u = ui(ri,0i) J ru 2 {r2,0 2 ) , 



(6) 



where the Ui(ri,0i) are the single- inclusion deformation 

fields implied by Eq. (j3j), with the bipolar coordinate 

1/2 

transformations T2 = (d 2 + r\ + 2dr\ cos#i) , cos 62 = 
(d + T\ cos^i) /r2, and sin 2 = (n sin#i) /r2, in which 
d is the center-to-center distance between the two inclu- 
sions. Note that, if the single-inclusion solution in Eq. ([3]) 
is considered up to some order n = TV Eq. ([6]) contains 
4 (2 N + 1) independent constants which must be fixed 
through the boundary conditions in Eqs. (j4]) and (j5j). 
Employing Eq. (|2j) , the elastic energy in Eq. (pQ) can then 
be written as a sum over the boundary terms 



Gi = -5 
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du 



d 



du \ 



dOi K h — V z u - K h u—V z u + ™— 



(7) 

at ri = Ri, where we have neglected terms independent 
of u. This expression is evaluated at each order in the 
Fourier-Bessel series in Eq. (j6j) using the orthogonality 
properties of trigonometric functions. Thus, the elastic 
energy in Eq. (pQ) reduces to an algebraic expression for 
the coefficients Af n and Bf n which, in turn, are pre- 
scribed by the boundary conditions in Eqs. ([!]) and (|5]). 

We determine the coefficients Af n and Bf n in Eq. (J6j) 
up to arbitrary N by expanding U2 in the membrane 
region surrounding inclusion 1 in terms of T\jd and im- 
posing on u the boundary conditions at T\ = R\ [13, Hl| , 
with a similar procedure for inclusion 2. Thus, finding the 
general solution in Eq. © which respects the boundary 
conditions in Eqs. (|4j) and (j5j) reduces to solving a sys- 
tem of 4(27V + 1) linear equations, which we achieve using 
standard methods [g9|. Substitution of the resulting ex- 
pressions of the coefficients into Eq. ([7]) yields, for a given 
TV, the interaction potential between two membrane in- 
clusions. The validity of this finite-order series solution 
is based on the assumption that very rapid angular varia- 
tions at large n in the membrane deformation field can be 



neglected, which we confirm for a given problem by sys- 
tematically including higher-order terms. We have ver- 
ified [70] this methodology for the analytic calculation 
of elastic interaction potentials in the case of cylindrical 
membrane inclusions of constant hydrophobic thickness, 
for which elastic interaction potentials were calculated 
previously 0, using computer simulations. 

We now apply the above methodology to estimate elas- 
tic interactions between MscL proteins. In the most basic 
scenario, MscL can exist in one of two states — open or 
closed — with the competition between these two states 
governed by membrane tension 43 -4Zj. A simple coarse- 
grained model of the cross-sectional shape of MscL is 



Ci(0i) = Ri [1 + €i COSS (0i - UJi)} 



(8) 



where €i parameterizes the magnitude of the deviation of 
the protein cross section from the circle, s denotes the 
order of the protein symmetry (ohgomeric state), and uoi 
captures the orientation of the protein [see insets in Fig.Q] 
for illustrations of Eq. JHJ)]. The physiologically relevant 
ohgomeric states of MscL remain a matter of debate 
l50L with tetrameric [ill , pentameric (5l| , and hexameric 
[721 ] MscL having been reported. Our methodology is able 
to handle all of these cases, but here we focus on pen- 
tameric MscL, which correspond to s = 5 in Eq. (Et with 
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Ri and ^ estimated from structural models 
Since these structural models of MscL suggest q <C 1, 
we concentrate on the weak perturbation limit of Eq. (|8j) 
and only consider the leading-order terms in break- 
ing rotational symmetry. Expansion of u at ri = Ci(0i) 
in ei [H, yields boundary conditions of the form in 
Eqs. ([!]) and ([5j). For simplicity we first consider the case 
of a constant Ui along Q(^), with a value of a so that 
Ui is of the same sign in the open and closed states of 
MscL, and then turn to the complementary case of direc- 
tional interactions induced by a varying Ui p3,B3. We 
choose [ill the parameter values for bilayer-MscL inter- 
actions as discussed in Refs. jHiMHHI and, in particular, 
set U!(0i) = along the bilayer-MscL interface. 

Figure [TJa) shows the elastic interaction potentials 
between pentameric MscL proteins in the open state 



4lL l52l - l55j . with similar curves for the corresponding 
closed state structure (5l|. Irrespective of the relative 
protein orientation considered, we find an energy bar- 
rier to the dimerization of two open (or closed) MscL 
proteins. For most protein separations, the tip-on ori- 
entation is energetically most favorable with the face-on 
orientation being least favorable. For very small d, how- 
ever, the face-on orientation becomes favorable over the 
tip-on orientation [see inset in Fig. DJa)]. Figure \Hb) 
demonstrates that elastic interactions between open and 
closed MscL proteins can also yield an energy barrier to 
dimerization. But, in contrast to Fig. DJa), MscL pro- 
teins now repel each other at very small d, and we obtain 
a pronounced minimum in the interaction energy at some 
optimal value of d which depends on the relative protein 
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FIG. 1: (color online). Elastic interaction potentials ob- 
tained from Eq. (0) for (a) open and (b) open and closed 
MscL. Bilayer-MscL interactions were parameterized [70[] as 
in Ref. [65[ but for a PC20 lipid bilayer with r = 0. Thick 
curves denote the three MscL orientations in the insets, while 
thin curves correspond to intermediate MscL orientations ro- 
tated by 7r/20 or 7r/10 as indicated by arrows. The molecular 
models of MscL in the insets are taken from Ref. |4l[ and the 
superimposed boundary curves were obtained from Eq. ((8]) 
with Ri & 3.49 nm and Ri « 2.27 nm [70] for open and 
closed MscL. 



orientation. Similarly as in the case of two open MscL 
proteins, the tip-on orientation is energetically most fa- 
vorable, and the face-on orientation least favorable, for 
most protein separations in Fig.QJb), with the tip-on ori- 
entation becoming least favorable for very small d. For 
all scenarios in Fig. [1] we find that G changes smoothly 
upon rotation from the face-on to the tip-on orientation. 

The analytic solution in Eq. (j6j) in terms of single- 
inclusion deformation fields suggests a simple qualitative 
explanation for the non-monotonic behavior of the MscL 
interaction potentials in Fig. [TJ In general, the single- 
inclusion thickness deformation field in Eq. (|3]) over- 
shoots considerably 
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26| at its first extremum when 
relaxing away from the inclusion boundary. Hence, if two 
channels induce thickness deformations of the same sign, 
there is a regime of intermediate d for which membrane 
deformations are amplified, yielding the energy barriers 
in Fig. [TJ At small the membrane deformation fields 
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FIG. 2: (color online). Elastic interaction potentials obtained 
from Eq. J7} for cylindrical membrane inclusions of radius 
Ri = 3.5 nm with the hydrophobic mismatch in Eq. ©. The 
lipid bilayer was parameterized [73] as in Ref. [SI] with r = 0. 
Thick curves denote the three inclusion orientations in the in- 
sets, while thin curves correspond to intermediate orientations 
rotated by 7r/45 as indicated by arrows. 



already overlap before reaching the first extremum, thus 
reducing the overall deformation footprint of the two 
channels and making dimerization favorable as implied 
by Fig. [TJ For membrane channels of distinct hydropho- 
bic thickness, however, moving the channels even closer 
together yields an additional regime in which the mem- 
brane has to deform strongly between the two channels, 
leading to repulsion between open and closed MscL at 
very small d as in Fig. QJb). In each of these regimes, 
the competition between tip-on and face-on orientations 
is governed by a complex interplay between the local 
strengths and associated membrane areas of overlapping 
deformation fields. 

The qualitative picture of elastic interaction potentials 
developed above for the boundary curves in Eq. (JHj) also 
applies to the complementary case of directional inter- 
actions induced by a varying hydrophobic thickness of 
membrane proteins 
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JJJ . To illustrate this scenario we 
consider cylindrical membrane inclusions with hydropho- 
bic mismatch 



Ui(0 i ) = U? + 5 i coBs(0 i -uj i ) 



(9) 



where U® is the average hydrophobic mismatch and Si 
is the magnitude of mismatch modulations [see insets in 
Fig. [2] for illustrations of Eq. flS)]. Figure [2] shows the 
interaction potentials associated with Eq. (j9j). A striking 
feature of Fig. [2] is that, depending on the relative inclu- 
sion orientation, membrane-mediated interactions can be 
attractive or repulsive at small d. In particular, if the 
periodic modulations in Eq. ([9]) are in phase at the clos- 
est approach of the two inclusions there is, again due to 
the oscillatory nature of the single-inclusion deformation 
fields, an energy barrier to dimerization as in Fig. HJa), 
with attraction at small d. Conversely, in the case of ori- 
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FIG. 3: (color online). Gating curves for a pair of MscL 
proteins to transition from the open-closed to the open-open 
configuration [see Figs. Hfb) and HJa)] for the lipid bilayer 
parameter values in Ref. [Sj. The gating curves at d = 
9 nm and d = 10 nm for different protein orientations are 
labelled as in Fig. [TJa). The molecular models in the insets 
are taken from Ref. |4l| , and we use the same boundary curves 
for open and closed states of MscL as in Fig. [T] 



entations in Fig. [2] for which the periodic modulations in 
Eq. (J9j) are out of phase at the closest approach of the two 
inclusions, opposing inclusion boundaries induce distinct 
thickness deformations which, similarly as in Fig. QJb), 
can lead to repulsion at small d. The interaction poten- 
tials in Fig. [2] change smoothly upon rotation from out- 
of-phase to in-phase orientations, and exhibit a minimum 
for intermediate orientations. 

Interaction potentials such as those in Figs. Q] and [2] 
can induce directionality in the cooperative function of 
membrane proteins. In particular, the gating character- 
istics of MscL with varying membrane tension are cap- 
tured 43 -45[ by the channel opening probability P Q — 
j(ag-tAA)] where p = i/k B T, in which k B 



V [i- 

is Boltzmann's constant and T is the temperature, and 
AG and AA are the free energy and area difference be- 
tween open and closed channel states. We only consider 
here contributions to AG due to thickness deformations 
of the bilayer membrane in Eq. flTJ), which were found 
previously [24], Ir32l - [65| to yield the basic phenomenology 
of MscL gating. Figure [3] shows P Q for pentameric MscL 
[4lL l47H57j with a neighboring pentameric MscL protein 
in the open state. We find that, compared to the case of 
non- inter acting MscL, membrane- mediated interactions 
can, depending on the protein separation and orienta- 
tion, shift the gating tension to higher as well as lower 
values. The magnitude of the effect of directional inter- 
actions in Fig. [3] is comparable to previously measured 
25l IHij shifts in gating tension due to modification of 
the membrane composition. As expected from the inter- 
action potentials in Fig. [TJ the cooperative gating tension 
associated with the tip-on orientation is lower than the 
cooperative gating tension associated with the face-on 



orientation for the protein separations in Fig. [3j with a 
smooth interpolation in P Q between these two limiting 
orientations. 

In summary, we have developed an analytic methodol- 
ogy for estimating directional interactions between mem- 
brane proteins for arbitrary protein separations. This 
methodology provides a bridge connecting the hydropho- 
bic shape of membrane proteins to the cooperative func- 
tion of membrane proteins induced by elastic interac- 
tions. Our methodology is general and can be applied 
to any membrane protein for which basic structural in- 
formation is available. On the basis of a simple model of 
the cross-sectional shape of MscL suggested by structural 
studies 41, 47 -5?}, we predict that bilayer- mediated in- 
teractions between MscL proteins yield a characteristic 
sequence of preferred orientations, and distinctive gat- 
ing tensions, upon dimerization of MscL. A combination 
of quantitative experiments on the spatial arrangement 
[24l [25| and gating tension [13, [H3 - [57| of MscL would 
be able to put these predictions to a direct experimental 
test. The approach developed here represents a first step 
towards a physical theory of how elastic interactions af- 
fect the molecular structure, organization, and biological 
function of proteins in the crowded membrane environ- 
ment provided by living cells. 
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